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Abstract 

In models with a low quantum gravity scale, one might expect sizable effects 
from nonrenormalizable interactions that violate the global symmetries of the 
standard model. While some mechanism must be invoked in such theories to 
suppress higher-dimension operators that contribute to proton decay, opera- 
tors that change baryon number by two units are less dangerous and may be 
present at phenomenologically interesting levels. Here we focus on AB = 2 
operators that also change strangeness. We demonstrate how to compute ex- 
plicitly a typical nucleon-nucleon decay amplitude, assuming a nonvanishing 
six-quark cluster probability and MIT bag model wave functions. We then 
use our results to estimate the rate for other possible modes. We find that 
such baryon-number-violating decays may be experimentally accessible if the 
operators in question are present and the Planck scale is less than ~ 400 TeV. 
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I. INTRODUCTION 



One of the exciting implications of extra spacetime dimensions with large radii of com- 
pactification is the possibility that the scale of quantum gravity may be brought down to 
TeV energies |IJ. While much effort has focused on understanding the experimental signals 
of the extra dimensions themselves (for example, through the effects of Kaluza-Klein ex- 
citations in the effective four- dimensional theory), much less has been said on the physics 
that originates at the cut off 0. This is understandable for two reasons. First, any detailed 
knowledge of Planck-suppressed operators in a theory with a low quantum gravity scale 
requires a complete theory of quantum gravity, which is not yet at hand. Alternatively, a 
general effective field theory approach, in which one includes all operators consistent with 
the gauge symmetries of the standard model and suppressed by powers of the cut off, leads 
to baryon- and lepton-number- violating effects far in excess of the experimental bounds [|3],[|] . 
The most common approach to this dilemma is to assume that some mechanism forbids the 
undesirable operators, and then to ignore the issue altogether. Here we will explore the pos- 
sibility that the mechanism responsible for maintaining proton stability does not forbid the 
complete set of operators that violate the global symmetries of the standard model. Assum- 
ing AB = 1 interactions are absent, higher- dimension operators that violate baryon number 
B by two units are far less problematic, and may be present in low Planck scale scenarios 
at a phenomenologically interesting level. Moreover, such operators are not suppressed by 
separating quarks and leptons in an extra dimension, as has been suggested as a remedy to 
the proton decay problem 0. In this letter we study the two-body double nucleon decays 
that are induced by AB = 2 operators and determine the sensitivity of existing experiments 
to the scale of the new physics. 

The idea that the most dangerous baryon-number- violating operators may be absent, 
while others are present is not at all a radical one. Consider any model in which baryon- 
number is promoted to a gauge symmetry and then spontaneously broken: higher-dimension 
operators that violate baryon number by AB units are induced in the low-energy theory, 
but the smallest possible value of AB is controlled by the charge of the Higgs field that is 
responsible for the spontaneous symmetry breaking [Q. From a low-energy perspective, the 
original continuous gauge symmetry is irrelevant, and one concludes that a discrete remnant 
is sufficient to eliminate the unwanted interactions. Such "discrete gauge symmetries" are 
known to be preserved by quantum gravitational effects 0, are well defined as fundamental 
symmetries 0, and arise in string theory ||. It is not hard to imagine scenarios in which 
operators that contribute to nucleon decay are forbidden by some residual discrete symmetry 
below the string scale, while operators of higher-dimension that violate baryon number 
remain in the low-energy theory M. 

Past interest in effective AB = 2 interactions has appeared in the context of grand 
unified theories [ID], R-parity- violating supersymmetry []TTJ| , and Planck-scale physics 



The relevant dimension-nine operators have been cataloged in the literature fll||T4|| , but have 
not been studied in their entirety. This is due in part to the difficulty in evaluating hadronic 
matrix elements, and the relatively large number of operators involved. The fact that 
each operator has an undetermined coefficient of 0(1) leads to an unavoidable theoretical 
uncertainty, and makes the value added in undertaking a complete analysis somewhat small. 
We will proceed by selecting a typical operator and decay process that is convenient for an 
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explicit matrix element evaluation (in fact, one that has never appeared in the literature); 
we then use this result to estimate the size of other accessible modes. 

Our canonical decay calculation will be for the process D — > K*K, where D represents 
a deuteron. We choose this mode for a number of reasons: (i) The underlying operator 
changes strangeness and is unconstrained by neutron-antineutron oscillation bounds. (We 
will have more to say about the relationship between the operator we consider and others 
that do not change strangeness in the final section.) (ii) the operator we consider contributes 
to the matrix element of this process via precisely one Feynman diagram, and (iii) the spin- 
flavor-color-spatial wave function of the initial state is easily cross-checked with studies of 
the deuteron structure that appear in the literature |0J. We compute the effective lifetime 
for this decay, and then extrapolate to other NN modes of interest. We will show that such 
two-body decays may be accessible if the Planck scale is less than ~ 400 TeV, and we discuss 
the dependence of our result on the flavor structure of the theory. 

How then can a significant number of quarks within a deuteron be annihilated by a 
contact interaction, if the quarks are spatially separated within the neutron and proton? 



II. DEUTERONOMY 

A favorable answer to this question is that the deuteron is not entirely made of a proton 



and neutron, but also includes a significant admixture of a six-quark (6q) cluster state |L5 
In the case at hand, the 6q cluster is a six quark state that is totally antisymmetric in color, 
spin, and flavor, with isospin and spin 1. Such a state is not factorable into a simple 
product of two color- singlet, three-quark states. 

When the constituents of the nucleus are far apart, a description in terms of neutrons 
and protons is accurate. The question is what happens to the material inside the nucleus 
when the pieces come close to each other. The quark cluster viewpoint is that if two nucleons 
come sufficiently close together, the quarks within them reorganize into a new state where 
each quark is in the lowest energy spatial state, and the color-spin-flavor part of the quark 
wave function is uniquely fixed by the requirement that it be totally antisymmetric, 
colorless, and of the desired spin and isospin. 

Quark clustering gained impetus [16] in explaining 3 He data in kinematic regions inac- 
cessible to scattering off a stationary nucleon and where contributions due to Fermi motions 
calculated in standard models were insufficient. Further, differences between quark distri- 



butions in 6q clusters and in nucleons provide one straightforward explanation [M of the 
nuclear EMC effect. 

The deuteron state may be expressed as a linear combination 

\D) = (1 - f^D, NN) + f^ 2 \D, 6q) , (2.1) 

where / is the 6q cluster probability. The crucial point is that there is no reason to expect the 
overlap of the quark spatial wave functions to be negligibly small in a 6q cluster. Providing 
/ is also non- negligible, we avoid the possibility of obtaining an uninteresting result due to 
wave function suppression. 



Values for / have been estimated from 6q cluster models of the nuclear EMC effect [17 



from calculations of the probability for nucleons to overlap using realistic deuteron wave 
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functions [18|; from descriptions of high energy SLAC electron-deuteron data at x > 1.0 [19|; 



from studies of the deuteron electromagnetic structure functions [M; and from calculations 



of fast backward tagged nucleons coming from deuteron breakup [EIJ. The estimates range 



between 0.01 and 0.07. We shall quote results assuming / = 0.01 for free deuterons, and 
scaling / appropriately for 6q clusters within the nuclear medium. 

It will be convenient for us to represent the spin-flavor-color structure of the deuteron 
6q state in a field-theoretic form. Letting represent an (anticommuting) quark creation 
operator, the highest projection spin state may be written 

\D,6q,S z = 1) = Ne a ^e 5< e ab e de e c} e ij e kl 

X <lLiilbfi\cAdkieiif^) > ( 2 - 2 ) 

where the Greek indices represent SU(3) color, a . . . e SU(2) flavor, and i . . . I SU(2) spin, 
and is a creation operator for a quark state satisfying 

{laai, Q/3bj} = ^a^ab^ij ■ (2.3) 

One may see by inspection that the state has the quantum numbers of the deuteron. Note 
that the SU(2) spin space corresponds exactly to the MIT bag model wave functions de- 
scribed below (so the reader should not think that our calculation is nonrelativistic). In this 
representation, it is straightforward (though tedious) to compute the normalization factor 
N- we find N = 1/(48^). 

III. MATRIX ELEMENT 

With the state defined, we now turn to one possible operator of interest 

O = [t£°C«£] [4 7 C4] \sl P Csl] e a , p e, Sr , (3.1) 

which contributes to the decay D — > K*K. Here, C represents the charge conjugation ma- 
trix, and all the fermion fields are right-handed, so that O is manifestly SU(3)xSU(2)xU(l) 
invariant. Since strange quarks are present only in the final state, the four remaining fields 
annihilate quarks in the deuteron, leaving two spectators. We represent the quark cluster 
component of the deuteron as well as the outgoing kaons as MIT bag model states (described 
in more detail below), and define the spatial origin as the point at which our AB = 2 oper- 
ator acts. When quark fields are replaced by ground state bag wave functions multiplied by 
appropriate creation or annihilation operators, the spin-flavor-color (SFC) matrix element 
may be factored from the spatial one. It may be determined by allowing 

6 = e Q ^e^ T e l h kl e mn u ai u^ k d 5l sl m 4 n (3.2) 

to act on the state Eq. Q2.2p , where the symbols now represent creation or annihilation 
operators rather than fields; one then takes the overlap with a similarly constructed two kaon 
state. We have computed the SFC matrix element by hand and by symbolic mathematics 
code, and obtain 

Vsfc = (K°;K*+,S Z = l\0\D,6q,S z = 1) SFC = -iVE . (3.3) 
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To evaluate the spatial part of the matrix element, we must take into account that the 
desired external states are eigenstates of momentum. To relate these to bag states, we 



suppose that the momentum eigenstates form a complete set, so that [22 

and similarly for the other states. The state with the subscript "B" is a bag state centered 
at R, spin variables are tacit, and the momentum eigenstate is normalized by 

(D(P')\D(P)) = 2E(2n) 3 6 3 (P - P') . (3.5) 

By inversion, 

2E r -> 

\D(P)) = — r Jd*R e- ip - R \D{R)) B . (3.6) 

The normalization condition determines 4>(P) = 4>{P), 

2 (P) =2E J d 3 re- lp ^I n (r) = 2EI n (P) , (3.7) 

where I n is the wave function overlap of two n-quark states centered at different points, 

I n (r)= B (D{-\f^\D{\r)) B . (3.8) 

For applications where the wave function is written as the product of n independent quark 
wave functions, one has 

Ur) = (h(r)T , (3-9) 

with 

h(r)= M-^^f)) b ■ (3-10) 
In the bag model, this may be written explicitly as 



h(r) = (4tt) / dz / pdp{u + u_ + l+l_R + ■ R_) (3.11) 

JO JO 

where, 

R ± = [(r/2Tz? + P 2 ] l/2 , R + -R- = (R 2 + + Rl-r 2 )/(2R + R_) (3.12) 

and where the factors Ui = u(Ri) and U = l(Ri) are the upper and lower components of bag 
wave functions 

m=(.; { ;l\ )=-^f..f* , . (3,3, 



Here x is a Pauli spinor, uo = x/R B = 2.0A3/R B , and N 2 R 3 B = x/[2(x - l)jo(x)] ~ 5.15. 

To calculate the matrix element for the decay D(P = 0) — > K*(p)K(—p), we apply the 
foregoing formalism to the quark cluster part of the deuteron as well as to the kaons. For 
simplicity, we assume the same bag radius R B ~ 1 fm for all states. We obtain 



M = (K(-p); K*{p), S z = 1\0\D(P = 0), 6q, S z = 1) 



2m D 2E 1 2E 2 



1/2 



x N 



(3.14) 



and 



N 



x 



B 



d 3 Rd 3 R 1 d 3 R 2 e ip{Zl - z ' 2) 

(K(Ri); K*(R 2 ),S Z = 1\0\D(R), 6q, S z 



1) 



B J 



(3.15) 



where E\ and E 2 are the energies of the K and K*, respectively, and the indirection is taken 
parallel to p. This may be re-expressed as 



N 



1 

8' 



TjsFC / d 3 Rd 3 R 1 d 3 R 2 e 



ip(Zi-Z 2 ) 



x /i(r 10 )/i(r 20 )K + l 2 Y(u lU2 + hhRt ■ R 2 ) 



(3.16) 



where u = u(R), I = l(R), and r i0 = \R± — R\. The factors of I\ come from the overlap of 
spectator quark wave functions. All the integrals in Eq. ( |3.14|) may be evaluated numerically. 
For pR B = 3.12, we find N = r] SFC (N 2 R 3 B /87v) 3 x 0.31, I 2 (p) = R 3 B x 0.41, and J 6 (0) = 
R 3 B x 0.13. 



IV. RESULTS 

Using the results from the previous section, it is straightforward to derive the decay 
width. In keeping with the proton decay literature, we will instead express our result in 
terms of a 'partial lifetime' (i.e. the lifetime if the branching fraction to the given mode 
were 100%). We find 

r(D -> K*K) = 2.18 x 10 7 yrs ■ (p^j / _1 (4.1) 

where M 5 is the dimensionful factor that suppresses the dimension- nine operator of interest. 
(While we identify M with the Planck scale in the present discussion, it is worth mentioning 
that our calculational framework is applicable to any model, e.g. R-parity violating super- 
symmetry, in which such effective operators are induced.) For a free deuteron (for example, 
in D 2 0), we set / = 0.01 and find r ~ 10 34 years for M ~ 293 TeV. For a deuteron within 
the O 16 nucleus, the quark cluster probability will scale as the nuclear density; for a reason- 
able estimate, / = 0.2, one finds r ~ 10 34 years for M ~ 395 TeV. Varying M by a factor 
of 2 changes the lifetime by 2 10 ~ 10 3 , which allows the result to vary from well below to 
significantly above the lifetimes usually associated with the maximum proton decay reach 
of Super-Kamiokande. 
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A more realistic estimate of the decay rate would also include flavor suppression factors 
that provide a small dimensionless number multiplying our operator. Without a model of 
flavor, however, we cannot exclude the extreme possibility in which Eq. ( |3.1| ) has no further 
suppression and all other AB = 2 operators are simply absent. On the other hand, we may 
consider how our results change if we impose a simple flavor ansatz: we could assume that 
fields of the first (second) generation are suppressed by a factor of A 3 (A 2 ), where A ~ 0.2 is of 
comparable size to the Cabibbo angle. (This could arise if a AB = 2 operator involving only 
third generation fields is allowed, and all others are generated from it via CKM-like rotations 
on the fields ||.) In this case, our previous estimate for r is extended by A -32 ~ 10 22 ; for 
decays in O 16 we now find r ~ 10 34 years for M = 2.3 TeV. With this ansatz, we can 
now also say something about the strangeness conserving operators that contribute to n-n 
oscillation. Estimates of the relevant matrix elements already exist fbl , from which we find 



- 10- yrs • (^) A- = 7.6 yrs. . (^) (4.2) 



which exceeds the best experimental bound, 3.8 yrs., 90% CL [23| for any value of M that 
could be plausibly identified with the Planck scale. 

Table 1 shows other possible nucleon-nucleon decay modes, and the multiplicative cor- 
rection that must be applied to our D — > KK* result to obtain an estimate for the lifetime. 
This factor takes into account the differing phase space, spatial wave function overlaps and 
flavor suppression following from our previous ansatz. We simply assume the SFC matrix 
elements are of comparable size. While we don't compute these matrix elements explicitly 
(given the far larger uncertainty from the proliferation of operator coefficients) the contri- 
bution of any given operator of interest may be evaluated explicitly using the approach we 
have presented. Observation of any of these modes in the absence of conventionally ex- 
pected nucleon decays would be a remarkable and unexpected sign of exotic physics not too 
far beyond directly accessible energies. 
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TABLES 



AS = 


(xA- 4 ) 


AS* = 1 


(xA" 2 ) 


AS = 2 


(xA°) 


7T7T 


2.54 


Kvr 


1.91 


KK* 


1.00 


piT 


1.40 


K*ir 


1.25 


KK 


1.49 


PP 


0.90 


Kp 


1.11 


K*K* 


1.10 


7T7] 


1.85 


K*p 


0.90 






PV 


1.06 


Kr\ 


1.40 






7T7]' 


1.18 


K*7] 


0.96 






prf 


0.96 


Krf 


0.96 






7]7]' 


0.94 


K*rf 


1.87 






rm 


1.32 











TABLE I. Lifetime correction factors for typical NN decay modes. The numbers do not include 
the flavor parameter A ^ 0.2, which should be multiplied in as indicated. 
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